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Solid-state diffusion in amorphous zirconolite 
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We discuss how structural disorder and amorphization affects solid-state diffusion, and consider zirconolite as 
a currently important case study. By performing extensive molecular dynamics simulations, we disentangle 
the effects of amorphization and density, and show that a profound increase of solid-state diffusion takes 
place as a result of amorphization. Importantly, this can take place at the same density as in the crystal, 
representing an interesting general insight regarding solid-state diffusion. We find that decreasing the density 
in the amorphous system increases pre-factors of diffusion constants, but not decreasing the activation energy. 
We also find that atomic species in zirconolite are affected differently by amorphization and density change. 
Our microscopic insights are relevant for understanding how solid-state diffusion changes due to disorder and 
for building predictive models of operation of materials to be used to encapsulate nuclear waste. 


I. INTRODUCTION 

A variety of radiation sources are created and used in 
science and technology. This includes an important area 
of energy generation in nuclear power stations, where ki¬ 
netic energy of fission products is converted into heat 
and electricity. In this and other applications, the en¬ 
ergy of emitted particles often has a two-fold effect. On 
one hand, this energy is converted into useful energy, 
by heating the material around the particle tracks. On 
the other hand, this energy damages the material and 
degrades its important properties, including mechanical, 
thermal, transport and so on, to the point that a ma¬ 
terial might lose its functional purpose. A problem in 
fission and future fusion reactor^, this issue is particu¬ 
larly acute in the process of safe immobilization of nu¬ 
clear waste, and constitutes one of the pressing issues 
that modern society face^l21. The issue is closely related 
to public acceptance of nuclear industry and therefore 
to the future of nuclear power. Regardless of the future 
of nuclear industry, the amount of accumulated nuclear 
waste is very large and is steadily growing while no ac¬ 
ceptable solution of its safe storage exists. 

Crystalline ceramics have been proposed to immobilize 
highly radioactive nuclear waste 2 6 . The main require¬ 
ment for the immobilization matrix, the waste form, is 
to prevent the radioactive isotopes from leaking and pol¬ 
luting the environment. This is perceived to be a very 
challenging requirement, given the high radioactivity of 
the nuclear waste and the long radioactivity time spans 
that extend from several thousand years for fission prod¬ 
ucts to several million years for actinides. As recently 
observecP, nuclear waste encapsulation involves “design 
problem the likes of which humanity had never before 
attempted, because it involved a time scale that required 
predictions of material and system behavior tens and 
hundreds of thousands of years into the future. Some 


perspective on the uniqueness of this temporal projection 
comes from the realization that the most ancient mon¬ 
uments of past engineering achievement, such as Stone¬ 
henge and the Pyramids, are barely five thousand years 
old.” With no direct testing possible over long periods of 
time, the decision of using a particular wasteform needs 
to be informed by indirect simulated experiments as well 
as detailed theoretical understanding of how irradiation 
affects the ability of the waste form to remain an effective 
immobilization barrier. 

Waste form alteration, the main concern for the safe 
encapsulation of nuclear waste, is a complex phenomenon 
involving diffusion, leaching and dissolution 7 ®. Under¬ 
standing the effects of irradiation on waste form durabil¬ 
ity is challenging because the pro cess is complex and in¬ 
cludes many mechanisms at work 2 3 7 l Generally, waste- 
form alteration as a result of irradiation can include 
chemical changes at the surface, reactions with water 
and other environmental agents, increase of solid-state 
diffusion of atoms in the bulk due to radiation damage, 
increased defect mobility, defects interaction and so on. 
The above processes can be system specific, yet have a 
common underlying mechanism, thermal diffusion. It 
is therefore important to understand diffusion in waste 
forms and its changes due to radiation damage. 

An important effect of irradiation on the waste form 
is the remarkable increase of diffusion as a result of 
radiation-induced amorphization 9 10 ! Seen in easily 
amorphizable materials, this effect has been thought to 
be either absent or reduced in materials resistant to 
amorphization by radiation damage 11 ! However, even 
most resistant materials such as ZrC >2 still show consid¬ 
erable damage in the form of a large number of well- 
separated point defect^. 

An interesting possibility is having a waste form which 
is amorphized by radiation damage yet still shows low 
levels of alteration and continues to be an effective im- 
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mobilization barrier. Zirconolite, ZrCaT^Oy, is one of 
the phases in SYNROC mixture of different ceramics 
each tailored to immobilize different ions present in the 
highly-radioactive nuclear waste' 1 ®, and is currently the 
preferred waste form by the UK National Decommission¬ 
ing Agency for immobilization of actinides. In more re¬ 
cent experimental 3 ^, zirconolite has been rendered X- 
ray amorphous by alpha-decay processes of Pu and, sur¬ 
prisingly and contrary to other materials, did not reveal 
significant chemical and physical alterations, witnessed 
by the absence of phase changes and microcracks even 
at fairly large volume increases. Furthermore, zirconolite 
maintain ed st rong elastic response and overall chemical 
durabilitjl™Sl, in contrast to other materials studied. 
At the same time, the aqueous durability of zirconolite 
is strongly affected by radiation damage^, with the evi¬ 
dence supporting diffusion-controlled ion exchange as the 
main mechanism of alteration of radiation-damaged zir¬ 
conolite. These results call for further detailed investiga¬ 
tion of the mechanisms involved in the alteration of this 
wasteform and diffusion mechanism in particular. 

According to current understanding, the increase of 
diffusion due to radiation-damaged system is due to the 
associated density decrease. Indeed, apart from rare ex¬ 
amples such as Si, radiation-induced structural changes 
and amorphization are accompanied by density decrease. 
The effect of density on the activation energy for diffu¬ 
sion U (the energy needed by an atom to jump from its 
surrounding cage to the neighbouring quasi-equilibrium 
location) has been well understood since the early work 
of Frenkel®. When interatomic separations are at their 
equilibrium values in a solid, U is too high for the jump 
event to take place during any feasible time. However, if 
the cage can increase its size (for example, due to ther¬ 
mal fluctuations) and open up a low-energy local diffu¬ 
sion pathway, the jump can proceed much faster. If A?’ 
is the increase of the cage size required for the jump to 
take place, U is equal to the work required to expand the 
change elastically, and is 

U = 8 nGrAr 2 (1) 

where G is shear modulus and r is the cage radius®. Note 
that the elastic energy to expand the sphere of radius r 
by amount Ar depends on shear modulus G only because 
no compression takes place at any point. Instead, the 
system expands by the amount equal to the increase of 
the sphere volume®, resulting in a pure shear deforma¬ 
tion. Indeed, the strain components u from an expanding 
sphere (noting that u —)-0 as r —> oo) are u rr = —2 b/r 3 , 
ugg = u 00 = b/r^ : giving pure shear uu = 0. 

From a theoretical standpoint, density is considered to 
be the main factor governing G and hence U®. Indeed, 
when density decreases as a result of radiation-induced 
structural changes, the interatomic interaction strength 
decreases, reducing G. Further, Ar decreases due to the 
increase of the cage volume. Therefore, if G is constant at 
constant density, Eq. (1) makes two predictions. First, 


U reduces due to density decrease, the widely anticipated 
result corroborated by more recent work on diffusion pro¬ 
cesses in glasses and viscous liquid d-®. Second, U does 
not change at constant density. 

In contrast to the density effect, the consequences of 
amorphization at constant density for diffusion are not 
understood. Indeed, if a structural change (e.g. amor¬ 
phization or large accumulation of point defects and their 
clusters as in ZrO^p- 1 ) takes place at constant density, the 
volume of the atomic cage around the diffusing atom does 
not change on average. However, the wide distribution 
of interatomic angles in the disordered structure leads to 
the appearance of both faster and slower local diffusion 
pathways (see Fig. 1), even if this structure is of the 
same density as the parent crystal, with the net effect of 
increasing the diffusion. This is an agreement with ex¬ 
perimental results reporting the decrease of G as a result 
of amorphization at the same density 19 . 

Hard to estimate theoretically, the combined effect of 
local fast and slow diffusion pathways at the same den¬ 
sity is important to understand from the waste form 
perspective 7 -l as well as from the general point of view 
of properties of disordered state. Indeed, since the early 
work®, the extent to which structural disorder affects 
system properties remains widely debated. For exam¬ 
ple, recent work aiming at resolving the long-standing 
debate about the origin of the Boson peak in the energy 
spectrum of glasses, has found that, contrary to earlier 
expectations, the difference of vibrational spectra and 
other important properties between the amorphous sys¬ 
tem and its parent crystal disappears once the densities 
of both systems are taken to be the sam^XI. Generally, re¬ 
cently accumulated evidence suggests that disorder leaves 
some important properties of the system surprisingly un¬ 
affected, but modifies other properties substantially 22 . 

In this paper, we address an important element of 
waste form alteration, the change of thermal diffusion due 
to radiation damage in the waste form. We use molecular 
dynamics (MD) simulations to study the change of diffu¬ 
sion of different atomic species in zirconolite as a result of 
structural disorder. Importantly, MD simulations enable 
us to disentangle the effects of amorphization and density 
increase on diffusion and discuss these effects separately. 
Such a separation is very hard to achieve in experimen¬ 
tal radiation-damaged samples. We find that a profound 
increase of solid-state diffusion takes place as a result of 
amorphization. Importantly, this can take place at the 
same density as in the crystal, representing an interest¬ 
ing general insight regarding solid-state diffusion. We 
find that increasing the volume in the amorphous system 
increases pre-factors of diffusion constants. We also find 
that atomic species in zirconolite are affected differently 
by amorphization and density change. Our findings are 
relevant for both understanding solid-state diffusion in 
the presence of disorder and for building predictive mod¬ 
els of operation of nuclear waste forms. 
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II. RESULTS 

There are several ways in which structural disorder 
can be introduced. Although direct simulation of colli¬ 
sion cascades gives detailed information about the nature 
of the damage, producing completely disordered struc¬ 
tures by multiple cascade overlaps is not practical, es¬ 
pecially for realistic high-energy events and large system 
sizes'. We prepared the amorphous structures by first 
melting the system at 5000 K, equilibrating the high- 
temperature liquid for 100 ps and subsequently quenched 
the liquid slowly to room temperature, 300 K. We note 
that amorphization by quenching the liquid can be dif¬ 
ferent from radiation-induced amorphization in several 
respect^ 10 l 24 l 25 ( however our main motivation is to ad¬ 
dress a fundamental question of how diffusion is affected 
by topological disorder in general. 

We use DL_POLY MD simulation packag^ 2 '^ and the 
system of 1056 atoms with the recent interatomic poten¬ 
tial fitted to zirconolite propertied^. We have simulated 
and relaxed three different zirconolite structures: crys¬ 
talline, high-density amorphous with density equal to the 
crystal and low-density amorphous zirconolite with 5% 
decreased density as in the e xper imental samples dam¬ 
aged by the radiation damagd^SJ. 10% of U atoms were 
introduced as a substitution for Zr atom, representing a 
typical waste load in waste forms. The interatomic po¬ 
tential for the U-0 interaction was taken from RefP^. 

Experiments on alteration of damaged waste forms 
are conducted at high temperature in order to ob¬ 
serve alteration and diffusion during laboratory time 
scale Pm iilMl, Similarly, we performed several MD 
simulations at temperature high enough to observe dif¬ 
fusion. Diffusion was observed as the linear time de¬ 
pendence of the mean square displacement, (r 2 ) = 6 Dt, 
where D is the diffusion coefficient. 

In Figure 2, we show representative (r 2 ) calculated for 
crystalline and amorphous zirconolite. We note that at 
short time, (r 2 ) crosses over from the oscillatory to the 
diffusive regime, acquiring the linear time dependence 
characteristic of diffusion (see Figures 2-3). We perform 
our subsequent analysis on the basis of the linear (r 2 ) a t 
diffusive regime at long times. We limit the analysis to 
5 ns in time since at longer times at high temperature 
we observe the signatures of recrystallization, witnessed 
by the appearance of peaks in Zr-Zr and Ti -Ti sublat¬ 
tices beyond the medium-range order. For the following 
analysis, we will use two equations for temperature de¬ 
pendence of diffusion coefficient D and hopping time r, 
the average time between two consecutive atomic jumps 
at one point in spac^® 


- u ) 
k B 77 

(2) 

U \ 

k B I7 

(3) 


where U is the activation barrier for the diffusion event, 
td is the shortest (Debye) vibration period of about 0.1 
ps and pre-factor D 0 is the high-temperature limit of the 
diffusion coefficient when r —> r B . 

Apart from O atoms, we observe no diffusion in the 
crystalline systems on the time scale of our simulations. 
On the other hand, we observe the diffusion of all atoms 
in the amorphous systems at the same temperature (see 
Figures 2-3). Importantly, this includes the diffusion in 
amorphous structures of the same density as the parent 
crystal. We therefore find that structural disorder at the 
same density increases the diffusion constant profoundly, 
an unexpected finding since the early work 16 , it was den¬ 
sity which was believed to govern the activation energy 
barrier for diffusion 18 . This represents our first main re¬ 
sult from this work. 

We next calculate main parameters of diffusion, U and 
Do, in amorphous systems and their change due to differ¬ 
ent density. We calculate U by fitting the data in Figure 
4 to Eq. (p| as InD = InDo — and show the results 
in Table l!""The calculated values of U in amorphous zir¬ 
conolite represent our next quantitative result enabling 
future prediction of how the waste form will operate dur¬ 
ing long time scales, as discussed below in more detail. 

In Table 1, we observe no appreciable differences of 
U between amorphous structures of different densities 
within the error due to scatter. The scatter is related to 
large fluctuations in the system at very high temperature 
we simulated, and is larger for the less numerous U atoms. 
We therefore find that the considered moderate density 
decrease, corresponding to t he ex perimental swelling of 
radiation-damaged zirconolitd^H, does not have a signif¬ 
icant effect on U. 

Since we do not observe diffusion in the crystalline zir¬ 
conolite apart from O atoms, we are unable to calculate 
U in the crystal and compare it to the amorphous sys¬ 
tem. However, we can estimate the lower bound of U in 
the crystal, Ui , using Eq. (j3j) : Ui = T h In ^, where r ; 
is the longest simulation time and Th is the highest tem¬ 
perature simulated in the crystal. Taking Ti = 2100 K 
and r; = 130 ns, we find Ui of about 3 eV. Ui is therefore 
of the same order of magnitude as typical U in crystals 
(typical U in crystals can be larger by up to about a 
factor of 2 for different atoms). 

A discernible trend in Figure 4 is the increase of pre¬ 
factor Dq at smaller density. We estimate this increase 
using two methods. First, we directly calculate the in¬ 
crease of Dq by fitting the MD data to In D = In Dq — v 
This gives in the range 2 — 3 for different atoms, 
where subscripts 1 and 2 refer to amorphous systems at 
larger (crystalline) and smaller density, respectively (see 
Table 1). Second, the range of can be estimated by 

subtracting In Di = In D 0 1 — ^ and In D 2 = In D 0 2 — ^. 
Using our earlier result that U are the same within the 
error, we find In S 23 = In S 3 . Then, the range of In 5 21 

by calculating In ^ at low and high temperature in Fig¬ 
ure 4, giving in the range similar to the first method, 


r = td exp 
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as follows from Table 1. We therefore find that the in¬ 
crease of diffusion pre-factors with system’s volume is 
appreciable. 

As mentioned above, O atoms stand out from the rest 
of atomic species in that their diffusion is seen in the 
crystalline zirconolite, representing an interesting hetero¬ 
geneity of atomic species in terms of diffusion. Weakly 
bound as compared to other atomic species, O atoms do 
not show discernible decrease of U as a result of amor- 
phization and density decrease (with the error quantified 
in Table 1 and present in Figures 2, 4). However, we ob¬ 
serve the increase of D 0 as a result of amorphization by 
the factor of about 3. Contrary to other atoms, however, 
this factor is not sensitive to density decrease (see Figure 
4 and Table 1). 


III. DISCUSSION AND SUMMARY 

There are several important insights from this work. 
First, we find that structural disorder and amorphiza¬ 
tion in particular, introduced to the system at the same 
density as the parent crystal, can result in a profound 
increase of solid-state diffusion. This is an interesting 
general insight in view that both earlier and current the¬ 
ories emphasize density as the main factor controlling the 
diffusion 16 18 . Our finding implies that the net effect of 
the appearance of slower and faster local diffusion path¬ 
ways in the disordered structure (see Figure 1) is diffusion 
increase. 

According to our result, faster local diffusion pathways 
dominate over slower ones in the same-density disordered 
structure. This is analogous to the well-known case of in¬ 
troducing equal amounts of harder and softer inclusions 
in an elastic matrix, with the result that the overall elas¬ 
tic response is mostly governed by the softer phase (in 
a simple example, this is illustrated by the inverse sum 
rule for bulk or shear moduli). The analogy is further 
relevant here because lower-density faster diffusion path¬ 
ways and higher-density slower pathways can be approxi¬ 
mately viewed as local regions that are softer and harder, 
respectively (see Figure 1). Then, the net effect is elastic 
softening of the system and hence smaller G , in agree¬ 
ment with experimental results 1 ®. This implies smaller 
U according to Eq. (1) and hence larger D according to 
Eq. (2), consistent with our current findings. 

We note here that contrary to the role of density which 
is well understood and quantified, the effect of disorder at 
constant-density is not amenable to a simple treatment 
or a model, but is important in a wider context of un¬ 
derstanding the essential differences between crystalline 
and amorphous systems. For example, it has been long 
thought that amorphous systems are notably different in 
terms of their vibrational spectra, yet recent evidence 
has found that this is not the case if the amorphous sys¬ 
tem is at the same densitj^. This is consistent with a 
wider picture emerging that many important thermody¬ 
namic properties of the system are insensitive to disorder 


due to the similarity of their spectra whereas other prop¬ 
erties, namely the transport properties such as thermal 
conductivity, are strongly affected^. 

Second, the increase of diffusion due to amorphization 
is important for the operation of waste forms. For ex¬ 
ample, under the typical waste load zirconolite will be¬ 
come amorphous from irradiation after about 1000 years. 
Being a very small fraction of time of operation of the 
waste form (100,000-1,000,000 years), this implies that 
solid-state diffusion will take place almost entirely in the 
amorphous state. This will take place with the associated 
increased diffusion constants that we have found in this 
work. Our results, and the calculated values of U in par¬ 
ticular, can therefore be used to predict the solid-state 
diffusion in the waste form during the most important 
period of its operation. 

Third, we have found specifically that the combined 
effect of amorphization and volume increase can affect 
both the activation energy and diffusion pre-factors, but 
that different atomic species can be differently affected. 
For example, the diffusion of most numerous and mobile 
O atoms is not affected by amorphization to the same 
extent as in other atomic species. 

In summary, by using MD simulations we have dis¬ 
entangled the effects of amorphization and density on 
solid-state diffusion, and showed that contrary to exist¬ 
ing theories, a profound increase of diffusion takes place 
as a result of amorphization at the same density. We have 
found that decreasing the density in the amorphous sys¬ 
tem increases the pre-factors of diffusion constants. We 
have also found that atomic species in zirconolite are af¬ 
fected differently by amorphization and density change. 
Our microscopic insights are relevant for understanding 
how solid-state diffusion changes due to disorder and for 
constructing predictive physics-based models aimed at 
predicting the performance of waste forms over long time 
scales. 
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FIG. 1. Schematic picture showing equivalent interatomic dis¬ 
tances and local diffusion pathways in the crystalline structure 
(left). The disordered structure with the same density as the 
parent crystal has a wide distribution of interatomic distances 
and angles and gives rise to both faster (solid line) and slower 
(dashed line) local diffusion pathways. 



FIG. 2. Mean-squared displacement of O atoms (averaged 
over all atoms) in crystalline and amorphous zirconolites of 
two different densities at 2000 K. 
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TABLE I. Activation energy U (eV) of atomic species in the crystalline and amorphous zirconolite at two different densities. 
Pre-factors Do (m 2 /s) are denoted as Doi (amorphous system at larger crystalline density), D 02 (amorphous system at smaller 
density) and -D crystal (pre-factor of difusion of O atoms in the crystal). The change of diffusion pre-factors is evaluated using two 
methods discussed in the text: calculating Do from direct fitting to calculated D (a) and estimating the range of Do at different 
temperature (b). Subscripts 1 and 2 refer to amorphous systems at larger (crystalline) and smaller density, respectively. 


Ca 

Zr 

Ti 

U 

O 

U (amorphous, same 2.80 ± 0.53 
density as crystal) 

3.05 ±0.73 

3.79 ±0.63 

3.20 ±0.86 

1.61 ±0.06 

U (amorphous, 2.80 ± 0.20 

smaller density) 

3.05 ± 0.71 

3.79 ±0.55 

3.20 ±0.54 

1.61 ±0.15 

U (crystal) 

(a) 




1.61 ±0.07 

ln^sa 0.91 ±0.09 

Q1 

in 

-'-''crystal 

(b) 

0.66 ±0.15 

0.76 ±0.13 

1.01 ±0.18 

0.0 ±0.2 

1.04 ±0.1 

In ^ 0.89 - 0.92 

In D ^ 01 

-'-'crystal 

0.62-0.71 

0.54- 1.00 

0.79-1.27 

0.94-1.11 



FIG. 3. Mean-squared displacement of Zr, Ti, Ca and U 
atoms (averaged over all atoms) in amorphous zirconolites 
of two different densities at 2000 K. 
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FIG. 4. ln(D) vs inverse temperature, ^ for different atomic 
species. The range of temperatures in the x-axis is 1400-2000 
K for O and 1870-2000 K for cation. The dashed lines are fits 
of both sets of points to the straight lines assuming that the 
slopes are the same within the errors present. 







